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Abstract 

We theoretically compute the interface thermal resistance between crossing single walled carbon 
nanotubes of various chiralities, using an atomistic Green's function approach with semi-empirical 
potentials. The results are then used to model the thermal conductivity of three dimensional 
nanotube pellets in vacuum. For an average nanotube length of 1 ^m, the model yields an upper 
bound for the thermal conductivity of densely compacted pellets, of the order of a few W/m- 
K. This is in striking contrast with the ultra-high thermal conductivity reported on individually 
suspended nanotubes. The results suggest that nanotube pellets might have an application as 
thermal insulators. 
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I. INTRODUCTION. 



Tremendous advancements in nanoscale heat transfer provide new alternatives to managing 
phonon transport inside crystalline nanosized devices. As the typical lenghtscale of systems 
reach those of the dominant phonon wavelengths, the transport of heat occurs in a coherent 
regime. Recently, various approaches have been developed to reconcile the wave behaviour 
of phonons with the diffusive and non-diffusive (ballistic) processes occuring inside confined 
solid state devicesi'^'^. Thermal management in short-scaled systems cannot be carried out 
by taking into account the classical Fourier law. Impurities, defects, Umklapp processes as 
well as roughness play a key role in phonon scattering: the flow of heat drastically differs 
from that at the macroscale^'^'^'^. 

For the past decade, thermal transport inside carbon nanotubes as well as inside semicon- 
ductor nanowires have been intensively investigated for their ballistic transport properties. 
Experimental mesurements on individual CNT as well as theoretical predictions confirmed 
a very high thermal conduct ivity^>^*ii'i^ii^. Nonetheless, even if single-walled CNT exhibit 
a ballistic regime in conducting heat, recent measurements showed that conductivity gets 
drastically reduced in composites^ii^ii^ii^. This is explained by the contribution of the inter- 
facial contact resistances. As the percolation treshold is reached by incorporating CNT, an 
interconnection pattern is formed and thus, each CNT gets connected to a network structure 
through different contacts depending on the mass concentration. These contacts are respon- 
sible for reducing the heat flux in the material. As the structure gets compacted, interface 
effects become predominant. 

A few works have reported on the impact of interface resistances from an atomistic 
approach^ii^i^iSiiSSiSSi;^. Ref. Q performed molecular dynamic simulations to predict the 
contact thermal conductance dependance due to the length of the overlap, and the spacing 
between two aligned overlaping CNT. This work showed that reducing the contact area in- 
creased the contact resistance. The predominance of the interface contact resistance over the 
nanowire/nanotube intrinsic resistance has also been reported from a theoretical point of view 
for the case of a planar substrate connected to a silicon nanowire at low temperatures2^>2^. 
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The physical mechanism of this very low conductance is explained in terms of a phonon 
density of states mismatch at the junction. 

Relevant numerical techniques used to retrieve thermal properties are mostly molecular dy- 
namics (MD), the Boltzmann transport equation, and atomistic Green's functions (AGF). 
MD is very well suited to solve high temperature problems because anharmonicities are 
usually included in the interatomic potentials, the many body interaction are by the way 
fully captured. When harmonic scattering (due to interfaces or defects for example) is pre- 
dominant, thermal transport is best described by atomistic Green's functions^!. Pertubative 
treatments can also be used to capture high temperature anharmonic phonon scattering 

In this paper, we first use the Green's function formulation of heat transport to obtain 
the phonon partial-transmission functions at the junction of two crossing overlaping CNT 
(CCNT) for different chiralities. Afterwards, we calculate the partial conductances at the 
junction and their dependence with CNT chirality. We then present a simple model involv- 
ing the previously calculated interface resistances, which predicts the thermal dependence 
of CNT composites conductivity as a function of the CNT mass density. We discuss the 
dependence of the thermal conductivity on the various physical parameters characterizing 
the pellet. An interesting result stemming from the calculation is that the thermal conduc- 
tivity has an upper bound, for dense pellets, which does not depend on the diameter of the 
nanotubes in the composite. This upper bound is several orders of magnitude smaller than 
the thermal conductivity of individual carbon nanotubes. We finally discuss the implications 
of these findings for novel technological applications. 



II. PHONON TRANSPORT AT THE JUNCTION BETWEEN CROSSING NAN- 
OTUBES. 



A typical crossing nanotube configuration is shown in Fig. [TJ The theory of harmonic 



interface thermal conductance calculation using Green's functions was developed in Ref. 
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It is straghtforward to generalize the theory in order to compute partial phonon transmissions 
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FIG. 1: Typical crossing configuration of two carbone-nanoubes in a nematic (compressed/aligned 
and non-isotropic) phase. 



between any two of the four nanotube edges. The case that we are interested in is the one 
where each nanotube is kept at a given temperature, equal on both edges. In such case, the 
thermal current across the junction is given by 



where fu{i) is the Bose-Einstein distribution corresponding to the temperature of the upper 
(lower) nanotube in the junction. The thermal conductance across this junction is obtained 
as (T = dJ/dT, and given by 





Ways to compute the transmission function have been explained in previous publications. 
Here we employed a three-region formula, given by 



where Fs = k34(g44 — g^)k43. Other ways of obtaining the transmission are also possible, 
like the two-region formula, for example, and the results are equivalent in all cases—. 

In order to compute the Green functions, a proper description of the system's atomic configu- 
ration and interatomic interactions is necessary. We have employed the Brenner potential to 
compute energies and forces between atoms belonging to the same nanotube^. Interactions 








(3) 



4 



between atoms in different nanotubes are not given by the Brenner potential. These are 
weaker interactions, and it is customary to describe them by a Lennard- Jones potential^. 
Several parametrizations can be found in the literature. We have opted for the one used 
by Zhong and Lukes-®, but results are not expected to change much if another reasonable 
parametrization were used. The parameters are e = AAlmeV and a = 0.228nm 

The stable atomic arrangements were obtained by energy minimization. A large supercell 
periodic array of perpendicularly crossing nanotubes was constructed for this purpose. The 
atoms in the last unit cells at the four edges of the two crossing nanotubes in the supercell 
were kept in a rigid structure corresponding to the relaxed isolated nanotube lattice. This 
is needed in order to seamlessly join the crossing part to the four semi-infinite contacts 
afterwards, for the Green's function calculation. The large size of the supercell ensures that 
relaxation is achieved within the crossing region. The two nanotubes' relative distance, 
crossing position, and angle around their axes, were varied until the lowest energy minimum 
was found. Afterwards, the interatomic force constants of the system were numerically 
computed by finite differences between the forces upon small (0.00001 A) displacements 
of each degree of freedom. Force constant matrices were independently calculated for the 
isolated nanotubes, to be used for the projection of the semi-infinite contacts. 

The self-energies corresponding to the four semi-infinite contacts were computed using a 
decimation technique'^^ . Then, the total Green's function was obtained from the previously 
calculated force constants for the crossing region, and the phonon transmissions were evalu- 
ated using Eqn. [3l 

We start by looking at the phonon transmission inside one individual CNT, in the presence 
of another CNT crossing it. Fig. [2] shows that the effect of the crossing is very small: the 
transmission function is very slightly reduced with respect to the integer transmission values 
for an isolated nanotube. In contrast, comparing Fig. [2] with Figs. [3] and H] shows that the 
phonon transmission between two different CNT is much lower than the one on the same 
CNT. Roughly speaking, this means that each nanotube can be considered to be at a constant 
temperature, and temperature differences exist only at the interfaces between nanotubes. 
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FIG. 2: Phonon Transmission functions through the same nanotube in the presence of the crossing 
one. 
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FIG. 3: Phonon Partial Transmission functions between two crossing nanotubes of identical chirality 
(Zigzag cases). 

Figs. [3] and |4] also reveal that at low frequencies, transmission is finite (non-zero.) This 
means that quantization effects will become observable at very low temperature, in the 
form of a linear temperature depencence of the thermal conductance. The interface thermal 
conductances evaluated from the transmission functions are shown in Fig. [5l as a function of 
temperature. Due to the strong suppression of the interface transmission above a few tens of 
THz, the interface thermal conductance temperature dependence is very weak above lOOK. 
This contrasts with the isolated nanotube thermal conductance, which increases noticeably 
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FIG. 4: Phonon Partial Transmission functions between two crossing nanotubes of identical chirality 
(Armchair cases). 

between lOOK and 300K^^. We will now use these computed thermal conductance values to 
estimate the thermal conductivity of CNT pellets, in the next section. 
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FIG. 5: Interface conductance for (5, 5), (6,6), (7,7), (8,8) - (5,0), (7,0), (8,0), (9,0) and (10,0) chiralities 
at 300K. Crossing configurations are set for two identical CCNT. 



III. THE THERMAL CONDUCTIVITY OF COMPACT CNT PELLETS 



A pellet is formed by a randomly arranged network of finite length nanotubes, intercrossing 
at various points along their lengths. As sketched in Fig. [6l there are three length parameters 
defining the pellet: / = nanotube length, d = nanotube diameter, and = density of 

junctions, with D ~ segment length between junctions. As we have seen in the section above. 
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placing a junction does not appreciably diminish the phonon transmission along the same 
nanotube. Thus, to a good approximation (even if the nanotubes are micrometers long, and 
have junctions every 10 nanometers) the temperature on each nanotube can be considered 
constant. For a macroscopic pellet, on average, the temperature of a nanotube whose center 
of mass is at Xj will correspond to Tj = Tcoid — VT ■ x,, where Tcoid is the temperature of the 
colder side. 

Partial Heat Flow 




FIG. 6: Length parameters defining the CNT pellet : I = nanotube length, d = nanotube diameter 
and D ~ segment length between junctions 



The heat flux is thus determined by the amount of heat traversing each individual junction 
across the pellet. Let us assume that the density of junctions per unit volume is l/D^, 
i.e. the average distance between junctions is roughly ~ D. The thermal current per unit 
cross section is then J/S = a{\ATjunct\)/D'^, where a is the junction's thermal conductance 
computed in the previous section, and {\ATjunct\) is the absolute temperature difference at 
the two sides of the junction, averaged throughout the pellet. 

In order to evaluate {\ATjunct\) , we average over all possible orientations Qi, Q2 of the two 
nanotubes, and contact positions Xi, X2 measured along the two nanotubes. Denoting the 
nanotube length by /, we have 



W/2 






dxidx2 






/:, 


d{cos 62) 




27r 




/- / dxi I dx2 I ^ I ^ o- 1^1 cos 6*1 — X2 cos 6^21 (5) 

/^^0.18/. (6) 



Combining this result with those in the previous paragraph yields the pellet's thermal con- 
ductivity, 

J 0.18/ 

Since the average distance between junctions cannot be smaller than the nanotube diameter, 
there exists an upper bound to the pellet's thermal conductivity upon compaction. In 
principle, this upper limit is larger the longer the nanotubes are. However, if the nanotubes 
are so long that they are effectively multiply bent in a noodle fashion, / in the equations 
is no longer the total nanotube length, but rather twice the main curvature radius of the 
nanotubes. It is possible to rewrite the previous model assuming that the nanotubes are 
closed circumpherences, yielding essentially the same results. 

An important fact that stems from the previous equation is that the thermal conductivity 
of the most highly compacted pellets should be roughly independent of the diameter of the 
nanotubes^^. To see this, it is enough to realize that the junction's thermal conductance 
is roughly proportional to the contact area, which in turn is proportional to the square of 
the nanotube diameter, d?. Since the minimum possible distance between junctions D is 
also proportional to the result in Eqn. [7] is independent of d. Obviously, this is only 
true at sufficiently high temperatures. For low temperatures the thermal conductance enters 
the quantized regime, and it no longer scales like the contact area (see discussion of Fig. [H] 
below.) 

It is convenient to express the pellet's thermal conductivity in terms of its density, so that 
results can be directly compared with future experiments. Realizing that there are two 
segments per junction, the pellet's density is 

P-2 , (8) 

where Pgraphene = 7.6 X 10^^ Kg/m^ is the surface mass density of graphene. The thermal 
conductivity is thus linearly proportional to the pellet's density, as 

0.18/ 
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The high compaction density hmit, assuming again that D ~ 2d, is pc 



Pgraphene 

2d 



To illustrate these results we have considered a typical nanotube length of l/xm. Fig. [7] 
shows the pellet's thermal conductivity as a function of its density, for different nanotube 
chiralities. A real composite will contain a range of chiralities and diameters, yielding an 
average of the individual results shown. Each of the lines shown is terminated at the highest 
possible density for nanotubes of that diameter. All the thermal conductivities in the dense 
limit are of the same order, and smaller than ~ 5 W/m-K. This very low upper bound for the 
thermal conductivity of the pellets contrasts with the high thermal conductivities (600 times 
larger) measured in isolated carbon nanotubes, which are comparable to that of graphite^*^. 




Mass Density [g/cm 



FIG. 7: Thermal conductivity as a function of CNT mass density for (5,5), (6,6), (7,7), (8,8), (5,0), 
(7,0), (8,0), (9,0) and (10,0) chiralities at 300K. Curves end-up when high density Umit is reached 
(horizontal Unes). 



Figure [8] shows the thermal conductivity in the high density limit, as a function of nanotube 
diameter, for various temperatures. Results at 300K and lOOOK are roughly independent 
of d, whereas at lOOK the results clearly decrease with increasing d. This is a result of 
thermal conductance quantization effects becoming important. In the low temperature limit, 
the thermal conductance across the junction is given by the zero frequency value of the 
transmission function times the quantum of thermal conductance. As shown in section II, 
the transmission function is always of the order of ~ 1 in the low frequency limit, regardless 
of diameter. This results in an interface thermal conductance that does not strongly depend 
on diameter at low temperature, which combined with Eqn. [7] results in the observed decrease 
as d increases. 
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FIG. 8: Thermal conductivity in the high density hmit as a function of the CNT diameters for 
ZigZag configuration (white symbol) and Armchair configuration (black symbol). Temperature is 
taken at lOK (circles), 300K (triangles) and lOOOK (squares) 



This suggests that, by carefully choosing the size of the nanotubes, it might be possible to 
tailor the high to low temperature thermal conductivity ratio of the pellets, with possible 
interesting applications in thermal switches and sensors. Conversely, it would be conceivable 
to use the measured temperature dependent thermal conductivity in order to determine the 
diameter of the pellet's constituent nanotubes. 

The very low thermal conductivity values of even the densest pellets suggest that there is 
room for carbon nanotube based materials in thermal insulation applications. Until recent 
times, many efforts have been devoted to trying to use nanotubes for enhancing materials' 
thermal conductivity. Given the inherently interesting structural and chemical properties of 
CNT's, low thermal conductivity carbon nanotube materials may be a still not well explored 
yet attractive subject of future research. 



IV. CONCLUSIONS 



We have computed the interface thermal conductance between crossing carbon nanotubes 
of various chiralities, using an atomistic Green's function formalism with interatomic force 
fields. We found that interfacial heat flow is mostly carried by low frequency phonons, while 
phonon transmission above a few tens of THz is strongly suppressed. As a result, the interface 
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thermal conductance depends weakly on temperature above lOOK. We have then developed 
a model of thermal conductivity for CNT pellets that incorporates the calculated interface 
thermal conductance values. Using this model, we have found that the thermal conductivity 
of densely compacted pellets is severly limited by an upper bound several orders of magnitude 
lower than the thermal conductivity of isolated CNT's or graphite. This suggests the possible 
applicability of CNT pellets in thermal insulation applications. 
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